Efficient vibrational state coupling in an optical tilted-washboard potential via 
multiple spatial translations and application to pulse echo 
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We measure the application of simple and compound pulses consisting of time-dependent spatial 
translations to coupling vibrational states of ultracold 85 Rb atoms in a far-detuned ID optical 
lattice. The lattice wells are so shallow as to support only two bound states, and we prepare 
the atoms in the ground state. The lattice is oriented vertically, leading to a tilted-washboard 
potential analogous to those encountered in condensed-matter systems. Experimentally, we find 
that a square pulse consisting of lattice displacements and a delay is more efficient than single-step 
and Gaussian pulses. This is described as an example of coherent control. It is striking that contrary 
to the intuition that soft pulses minimize loss, the Gaussian pulse is outperformed by the square 
pulse. Numerical calculations are in strong agreement with our experimental results and show the 
superiority of the square pulse to the single-step pulse for all lattice depths and to the Gaussian 
pulse for lattice depths greater than 7 lattice recoil energies. We also compare the effectiveness of 
these pulses for reviving oscillations of atoms in vibrational superposition states using the pulse-echo 
technique. We find that the square and Gaussian pulses result in higher echo amplitudes than the 
single-step pulse. These improved echo pulses allow us to probe coherence at longer times than in the 
past, measuring a plateau which has yet to be explained. In addition, we show numerically that the 
vibrational state coupling due to such lattice manipulations is more efficient in shallow lattices than 
in deep lattices. The coupling probability for an optimized single-step pulse approaches 1/e as the 
depth goes to infinity (harmonic-oscillator limit), while in shallow lattices with large anharmonicity, 
the coupling probability reaches a maximum value of 0.51 for a lattice depth of 5 recoil energies. 
For square and Gaussian pulses the coupling in the lattice is even stronger, reaching maxima of 0.64 
at 6 recoil energies and 0.67 at 5 recoil energies, respectively. 

PACS numbers: 03.67.-a, 32.80.Lg, 32.80.Qk, 85.25.Cp 



I. INTRODUCTION 

Cold atoms in optical lattices are of great interest for 
the study of numerous phenomena ranging, from predic- 
tions of condensed-matter theory P, 0,H, 0, @ to quan- 
tum chaos [1, 0, H, [H to quantum information protocols 
03 EH, EH- The real-time adjustability of the parameters 
of an optical lattice make it an attractive system for ex- 
ploring the controllability of quantum systems. Control 
of coherent quantum states is a topic of growing excite- 
ment, and plays for instance a crucial role in the field of 
quantum information processing. In practice this control 
needs to be exercised in real world systems where the in- 
teraction Hamiltonian, and particularly the interactions 
with the environment, may not be known precisely and 
where the control parameters are limited in number. We 
use cold atoms in a vertically oriented optical lattice as 
a prototype system and study how well coherence can 
be controlled in the resulting tilted-washboard potential, 
which is also relevant to a number of other physical sys- 
tems, when the control parameters are restricted to the 
(time-dependent) phase and amplitude of this potential. 
Our particular focus is on controlling and maintaining 
coherence between quantized vibrational states in such a 
potential. 

Quantized motion of atoms in optical molasses and lat- 
tices was first observed in the early 1990s [HI, EH EE EH 
using methods of high-resolution spectroscopy of reso- 



nance fluorescence, and stimulated Raman spectroscopy. 
Manipulations of the optical potential have been used 
to create coherence between vibrational states; examples 
are breathing-mode oscillations [ItJ through a sudden in- 
crease of potential depth and by parametric drive, Rabi 
type oscillations [18l ] through periodic phase modulations 
of the lattice beams, and collapse and revival of wave- 
packet oscillations [19| after a sudden shift of the lat- 
tice. The experiments where sudden shifts of the lattice 
were used to couple vibrational states [H, H3| were per- 
formed in deep lattice potentials (U a > 350-Er), where 
the coupling between the lowest two vibrational states 
approaches the harmonic limit. The periodic phase mod- 
ulations [HI in a very shallow lattice (U ~ 6Er) did 
not result in large coupling cither. In our experiments 
in moderately shallow lattices (18 < U < 20)Er we 
demonstrate how a combination of sudden shifts of the 
lattice and delays may result in large couplings. Non- 
adiabatic shifts of the potential have also been used in 
a pulse-echo experiment [2(| to revive wave packet os- 
cillations. In a pulse-echo experiment, one would ide- 
ally want to completely invert the state populations by 
applying a 7r-pulse for a perfect revival of oscillations. 
Unfortunately, the width of the bands as well as inhomo- 
geneous broadening have prevented any one from demon- 
strating ideal 7r-pulses for all lattice quasi-momentum 
components simultaneously. Echo pulses consisting of 
lattice displacements and delays presented in this paper, 
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while not perfect 7r-pulses, do result in partial inversion 
of state populations, and the strongest echo signals ever 
observed in such systems. There have also been propos- 
als [U, [12] for laser-driven vibrational state transitions 
by adiabatic passage along dressed energy levels for ions 
in harmonic and hard-wall potentials. Experimentally, 
Raman transitions have been successfully used to create 
Fock motional states for ions trapped in a harmonic po- 
tential 23] and for cold atoms in far-detuned ID optical 
lattices [2J]. In fact, a displacement of the lattice acts 
like a pair of Raman pulses though it does not result in a 
complete population inversion. Here, we present a study 
of the extent to which such lattice manipulations alone 
can result in inversion of populations between the lowest 
two bands. It is interesting to note that our tilted op- 
tical lattice at depths supporting two or three bands is 
very similar to anharmonic three-level potentials in su- 
perconducting qubits [H, H^, [27, 28]. In these systems, 
it is of fundamental interest to learn how manipulations 
of the potential parameters can be best used to control 
the quantum states of the system. 

In our experiment, first we study the efficiency of three 
different pulses consisting of lattice displacements and de- 
lays in coupling the atoms in the first band to the second 
band. Next we compare the efficiency of these pulses in 
reviving the oscillations of atoms that are in a superposi- 
tion of vibrational states through a pulse-echo technique. 
Throughout this paper we refer to the first experiment 
as the coupling experiment and to the second experiment 
as the echo experiment. 

In section II we start with a description of the experi- 
ment, then in section [II Al present our results in the cou- 
pling experiment and compare them with simulations. In 
section TlIB I we compare echo amplitudes after the three 
experimental optimized pulses. The details of our nu- 
merical calculations are provided in section III. 



II. EXPERIMENT 

We trap 85 Rb atoms in a lattice detuned 25GHz above 
the D2 transition line, F = 3^F 1 ' = 4(A = 780nm). 
For typical depths of 0.6 — 0.7^,K, this results in a pho- 
ton scattering rate of 16Hz, equivalent to a lifetime of 
w 60ms, which is much longer than the time scale of 
the experiment. The vertical lattice is formed by two 
laser beams with parallel polarizations intersecting at 
an angle of 9 = 49.6°, resulting in a lattice spacing of 
a = ^- = 0.93/ito, where = ^sin(|) is the lat- 
tice vector. We load the lattice from optical molasses 
with the final temperature of ~ 10/i-K" and atomic den- 
sity of « 10 10 atoms/cm 3 . The molasses is kept on in 
the presence of the lattice (depths of 18 — 20Er, where 
E R = h 2 k 2 L /2m = h(2n x Q8bHz) is the effective recoil 
energy for our lattice) and switched off quickly. Approxi- 
mately 10% of the atoms are loaded to the lattice, which 
corresponds to a density of « 5 x 10 5 per lattice plane; 
there is no transverse confinement, and this density is low 



enough that we may neglect interactions between atoms. 
Due to the short coherence length of atoms in optical mo- 
lasses (« 60nm at lOfiK ) there is no coherence between 
wells. 

Figure [2a) is a plot of periodic potential and energy 
bands for a lattice depth of U = 18E R and Fig. QJb) 
shows the same potential in presence of gravity, where the 
force of gravity adds a tilt to the lattice with a magnitude 
of 2.86-Er per lattice site. 

The atoms are prepared in the lowest vibrational band 
of the lattice by lowering the lattice depth adiabatically 
to the point that only one band is supported by the lat- 
tice. The lattice is then kept at this depth for 4ms during 
which the unbound atoms accelerate downward and es- 
cape the lattice region due to gravity. The lattice depth 
is adiabatically increased back to a value that supports 
two bands, which corresponds in our experiment to po- 
tential heights of U = 18 — 20Er. With these operations 
we are able to prepare > 95% of the atoms in the low- 
est band while the remaining < 5% are either lost from 
the lattice or are excited to the second band. The Bloch 
oscillation period in this lattice is « 530^ts, but we do 
not observe Bloch oscillations because of the incoherent 
filling of the lattice, which populates all quasi-momenta 
with the same probability. Due to the tilt of lattice, the 
third band at lattice depths of U — 18 — 20Er is not 
bound (see Fig.QJ. In practice, this means that the first 
two states remain bound throughout the duration of our 
experiment, while the third state escapes via Landau- 
Zener tunneling in w 830/is; see section Mil for more 
details. 

The depth of the lattice is calculated from the mea- 
sured oscillation frequency between the lowest two bands. 
The technique for this frequency measurement will be 
described later; sample oscillations for a lattice depth of 
U = 20Er are shown in Fig. [6j The measured average 
oscillation period in this figure is T ave = (187±2)/is and 




qx qx 

Figure 1: (a) Lattice potential U(x) = U sin 2 (khx) for the 
depth Uo = 18-Er and the first three bands as a function of 
dimensionless parameter qx = xir/a, where a is the lattice 
spacing, (b) Schematic of the tilted lattice showing U(x) + 
mgx and the quasi-bound states of the Wannier-Stark ladder. 
The energy levels are the averages of band energies and they 
differ from site to site by 2.86Er, where 2.86Er is the tilt due 
to gravity per lattice site. 



3 



(a) Experiment 



0.8 


- h 




* Single-step 
■ Square 

* Gaussian 


0.6 


• 


% 








0.4 










0.2 






i 

i * 






n 









0.2 0.3 0.4 0.5 

Ax/a 



1 


(b 


Model 


0.8 


MA ' 






* Single-step 
■ Square 

* Gaussian 


0.6 


■ 












0.4 




".' 

■ a * 

■■■ 








0.2 




■ 













■ 









0.1 0.2 0.3 0.4 0.5 
AX/O 



(a) Experiment 



0.4 






* Single-step 
■ Square 

* Gaussian 


0.3 












0.2 


i 

.• 


h 


¥ % 


i 




0.1 


I 4 


f 


% 


s 















mm 



(b) Model 



Single-step 

Square 

Gaussian 



0.1 0.2 0.3 0.4 0.5 OA 0.2 03 OA 0.5 



Ax/a 



Ax/a 



Figure 2: (Color online) (a) Experimental Pn as a function 
of displacement after application of single-step, square and 
Gaussian pulses for a lattice depth of Uo = 18Er. The tem- 
poral widths of the square and Gaussian pulses are the exper- 
imental optimized temporal widths, (b) Calculated Pn for 
the three pulses, with calculated optimized temporal widths 
for the square and Gaussian pulses. For experimental and 
calculated optimized pulse parameters see table [T] 

the rms width of the decaying Gaussian is (258 ± 6)^s. 
The oscillation measurement is an ensemble measure- 
ment and the measured frequency is the average of a 
distribution of frequencies due to the transverse Gaus- 
sian profile of the laser beams forming the lattice. From 
the measured lattice beam detuning and peak intensity 
we calculate the potential at the center of the lattice to be 
(24±2)Pr. This lattice depth is in the range of depths we 
calculate from the Gaussian rms width of the decaying 
oscillations. More details on the stages of preparation, 
operations and measurement including the adiabatic fil- 
tering and displacements of the lattice can be found in 
[29l ] . The data presented in this paper are for two lattice 
depths; in the coupling experiment the lattice depth is 
U = ^8Eji and in the echo experiment U = 20Pr. 

A. Coupling Experiment 

The three pulses used in this experiment are single- 
step, square and Gaussian profiles in time. The single- 
step pulse is a single lattice displacement. The square 
pulse consists of a displacement, delay and a displace- 
ment in the opposite direction. The Gaussian pulse con- 
sists of small displacements of the lattice in time steps 
of 5/is, which combined together approximate a Gaus- 
sian shaped pulse. The lattice is displaced by changing 
the relative phase of the laser beams forming the lattice 
through acousto-optic modulators (AOMs). The relevant 
parameters of each pulse are optimized by maximizing 
the number of atoms transferred from the first band to 
the second band: the single-step pulse is optimized with 
respect to the amplitude of single lattice displacement; 
the square pulse is optimized with respect to the ampli- 
tude of displacements and the delay between the displace- 
ments; and the Gaussian pulse is optimized with respect 
to the maximum amplitude of the lattice displacement 



Figure 3: (Color online) (a) Experimental P12 as a function 
of displacement after application of single-step, square and 
Gaussian pulses for a lattice depth of U a = 18Er. (b) Calcu- 
lated P12 for the three pulses. Experimental and calculated 
temporal pulse widths are the same as in Fig. [2] 

and the FWHM temporal width of the Gaussian. 

Starting with essentially all atoms in the lowest band 
(n = 1), we measure the number remaining in band 
1 and the number transferred to band 2, to infer cou- 
pling probabilities Pn (the probability of remaining in 
band 1) and P12 (the probability of being transferred to 
band 2). We compare these results with values calcu- 
lated according to Pi„ = (l/27ra) /^" a dg|ci n |g, where 
\ c in\1 = \(^n,q(x)\0\ 1 i'i tq (x)}\ 2 is the square of the ma- 
trix element of the operator O between Bloch states of the 
same quasi-momentum in the first band and the states 
in the higher bands. After preparing atoms in the lowest 
band we apply each of the three pulses and measure Pn 
and Pi 2 as a function of lattice displacement. The com- 
parisons are made for a lattice depth of U — 18Pr with 
the average oscillation period of T12 = 200/xs between 
the lowest two bands. Both the measurements and the 
calculations are done for lattice displacements in steps 
of 0.0139a, equivalent to a relative phase shift of 5° be- 
tween the lattice beams. The pulses in the experiment 
and the model have been optimized independent of each 
other. Figure [2Ja) compares the experimental Pn as a 
function of displacement due to the three pulses. In the 
case of square and Gaussian pulses the optimized exper- 
imental temporal pulse widths have been used (see ta- 
ble |T] for optimized pulse parameters), which match the 
temporal widths we find from theoretical optimization. 
In figure H^b) we show the calculated Pn as a function 
of displacement after the three pulses. There is a very 
good agreement between the experiment and the model 
forPn. 

Figures [3ja) and[3jb) show the experimental and cal- 
culated P12, the coupling to the second band after each 
pulse as a function of displacement. The temporal pulse 
widths are the same as in Fig. O In the case of the single- 
step pulse, we find the largest transfer to the second band 
when the lattice displacement is 0.22a = 205nm, where 
a = 0.93/im is the lattice spacing. For the square pulse, 
the optimum occurs for a delay of 0.35Pi2, and a dis- 
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Experiment 
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Optimized Pulse 


A pu l se 




P12 




\^pu I s e 


P12 


Single-step 


(0.22±0.04)a 




0.26 ±0.01 


0.250a 




0.391 


Square 


(0.14 ± 0.02)o 


(0.35 ± 0.03)Ti2 


0.33 ±0.01 


0.154a 


O.35OT12 


0.477 


Gaussian 


(0.16 ±0.02)o 


(0.29 ± 0.03)Ti2 


0.31 ±0.01 


0.186a 


O.294T12 


0.466 



Table I: Optimized pulse parameters (spatial amplitude A pu i ae and temporal width W pu ise) and population transfer to the 
second band P12 after each pulse for a lattice depth of U a = 18Pr, with all atoms initially in the lowest band. The temporal 
width reported for the Gaussian pulse is the FWHM width. A pu i ae is given in terms of the lattice spacing a = 0.93/im, and 
Wp U i se in terms of average oscillation period between the lowest two bands T12 = 200/^s. The error bars for ApuUe and W pu ise 
indicate the uncertainties in our determination of the optimal pulse parameters, while those for Pi 2 are the direct measurement 
errors. 



placement of 0.14a = 130nm. We find the optimum pa- 
rameters for the Gaussian pulse to be a maximum dis- 
placement of 0.16a = 149nm and a FWHM temporal 
width of 0.294Ti2. Table[T]is a summary and comparison 
of experimental results with our numerical model, which 
is described in detail in section ITTT1 The temporal widths 
of the experimental pulses match the model exactly, while 
the amplitudes of the optimal pulses in the experiment 
are smaller than the ones calculated in the model. We 
believe this is mainly due to the lattice inhomogeneity. 
Despite this difference, the optimized pulse parameters 
are remarkably similar for the experiment and the model. 
The main difference is in the maximum couplings to the 
second band, which in the case of experiment are « 70% 
of the calculated ones for the three pulses (see Fig. [3] 
and Table UJ). The loss from the lattice after each pulse 
due to couplings to higher bands not bound in the latice 
is larger in the case of experiment than in the calcula- 
tion, which we attribute to experimental imperfections. 
This difference becomes more pronounced as the magni- 
tude of lattice displacement increases to the maximum of 
Ax/a = 0.5 (see Fig.©. 

Overall, both the experimental results and the calcu- 
lations show the same trend in the performance of the 
pulses, i.e., both show that the square pulse results in 
higher coupling to the second band than the single-step 
and Gaussian pulses. The square pulse is an example of 



coherent control. The main feature of a coherent control 
scheme is the ability to steer a process toward a desired 
outcome. In the simplest case, two different processes 
can lead to the same final state and interference between 
these processes allows one to manipulate the branching 
ratio into this final state by controlling the relative phase 
of the two processes [30| . For the pulses discussed in 
this paper, both the square and Gaussian pulses intro- 
duce a phase factor in the superposition state through 
time delays. Consider the simpler of the two pulses, the 
square pulse, in which the time delay is introduced in 
one step. Each displacement in the square pulse couples 
states 1 and 2, analogous to a beam splitter in a Mach- 
Zehnder interferometer. The appropriate choice of the 
relative phase between the two displacements maximizes 
the transfer to state 2. Figure [4] shows the integrated cou- 
plings P12, Pn and loss vs. scaled delay after the square 
pulse at the fixed displacement of Ax = 0.16a. The plot 
clearly shows the coherent oscillations of the populations 
as a function of delay, demonstrating that the final pop- 
ulations are determined by interference between the two 
paths. 

Not only can phase be used to enhance transfer into 
the desired state via constructive interference, but it can 
simultaneously minimize loss through destructive inter- 
ference into the undesired state. As Fig. 0] shows, the 
coupling Pi 2 is peaked at t < 0.5, while the coupling to 
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Figure 4: (Color online) Calculated couplings P12, P11 and 
loss vs. scaled delay for the optimized square pulse at fixed 
Ax = 0.16a in a lattice depth of U a = 18Pr. The ground and 
excited states populations coherently oscillate as a function 
of delay. 



Figure 5: (Color online) Calculations: Loss vs. coupling Pi 2 
after the square, Gaussian and single-step pulses in a lattice 
depth of Ua = 18.Er. For the same coupling P12, the square 
pulse results in less loss than the single-step and Gaussian 
pulses. 
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higher levels (which result in loss) peaks later. This en- 
ables us to choose an optimal time delay for high coupling 
P12 but low loss. This result is analogous to techniques 
which have been used in other exp eriments to optimize 
coupling while minimizing loss [3l| . 

To elucidate the utility of such coherent control for 
optimizing the branching ratio, we plot in Fig. [5] loss 
versus coupling P12 for the single-step, square and Gaus- 
sian pulses. The displacement magnitude increases as the 
curves are traced counterclockwise from zero. For the 
same coupling P12, the square pulse results in less loss 
than the single-step and Gaussian pulses. As P12 grows 
to 0.39 (the maximum achievable with a single-step pulse 
in a lattice depth of 18Er), the loss due to the single- 
step pulse becomes nearly 3 times larger than that due 
to a square pulse; for larger P12, the single-step pulse is 
insufficient, regardless of loss. The square and Gaussian 
pulses have similar effect but still the square pulse out- 
performs the Gaussian pulse. The surprising result here 
is that the square pulse involving a pair of hard pulses 
does a better job of minimizing loss for a given transfer 
probability into the excited state than the Gaussian soft 
pulse despite the frequency selectivity of the Gaussian 
pulse [Hj]. 



B. Pulse-echo Experiment 

To measure the efficiency of each pulse in reviving os- 
cillations, we first prepare atoms in the lowest band, then 
displace the lattice by a/6 = 155nm, where a is the lat- 
tice spacing. The displacement of the lattice puts atoms 
in a superposition of the lowest two bands. This dis- 
placement magnitude is chosen as a trade-off between 
optimum coupling and loss of atoms from the lattice due 
to a single shift. (We have verified experimentally that an 
initial displacement of this size results in larger original 
oscillation amplitude and a larger echo amplitude than 
any other initial displacement.) We then map the phase 
evolution of the superposition state into the population 
evolution of the lowest band |1) by another displacement 
of the lattice. Mathematically, these operations can be 
written as 



P f (t) = (l\D\\a\)R(ujt)pR\u)t)D(\c 



ID 



(i) 



where p is the initial density matrix describing the su- 
perposition state, -D(H) = exp{— iP\a\/K) is the dis- 
placement operator, and R(u)t) = exp(—iHt/h) is the 
time evolution operator. This technique is analogous to 
the Ramsey technique of separated oscillatory fields [12] • 
The number of atoms left in the lowest band after the 
above operations is measured by lowering the lattice un- 
til all other bands become untrapped, and counting the 
remaining atoms using fluorescence imaging (29j . 

Figure [6] shows a typical measurement of the evolu- 
tion of Pf(t). The oscillations have an average period 
of T ave = (187 ± 2)ps with rms width of (258 ± 6)/xs. 
We believe much of the dephasing of oscillations is due 
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Figure 6: (Color online) Experimental data: Evolution of 
Pf(t) as a function of time. The curve is a fit to data (si- 
nusoid with Gaussian decay). T ave = (187 ± 2)ps and rms 
width of Gaussian is (258 ± 6)ps 



to the transverse Gaussian profile of the laser beams, 
which results in an inhomogeneous distribution of lattice 
depths and therefore an inhomogeneous distribution of 
oscillation frequencies. The oscillations can be restored 
by inverting the populations of the two states, i.e., ideally 
by applying a 7r-pulse. This results in a revival of oscilla- 
tions known as the echo signal. Spin-echo [33j in nuclear 
magnetic resonance and photon-echo [34| in optical res- 
onances are famous examples of this type of revival. 

The echo pulse is applied when there are no observable 
oscillations, usually after 3-4 periods of oscillations. For 
an echo pulse centered at time t , a revival of oscilla- 
tions or the echo signal centered at time 2t is observed. 
Figure a) shows the full trace of the lowest band popu- 
lation evolution after applying different pulses. All pulses 
are centered at t — 1040/is. The initial oscillations on 
each curve are due to the echo pulse itself since it is not a 
perfect 7r-pulse and creates its own coherence. These os- 
cillations dephase on the same time scale as the original 
oscillations due to the same inhomogeneous broadening. 
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Figure 7: (Color online) (a) Full echo curves after applying 
each pulse centered at time to = 1040/is. All echo signals 
are centered at 2to = 2080/is. The initial oscillations are due 
to the echo pulse itself, (b) Subtracted echo curves where 
the initial oscillations due to pulses measured after 4ms have 
been subtracted from the curves in (a). All curves have been 
displaced vertically for clarity. 
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Experiment 




Model 


Pulse 




Wp U l 3e 


Echo amplitude 


Pn>3 


Pn>3 


Single-step 


(0.29±0.05)a 




0.19 ±0.01 


71% 


52% 


Square 


(0.17±0.03)a 


(0.40 ± 0.05)Ti2 


0.23 ±0.01 


70% 


48% 


Gaussian 


(0.17±0.03)a 


(0.31 ± 0.06)Ti2 


0.24 ±0.01 


55% 


39% 



Table II: Echo amplitude and number of atoms lost due to each pulse centered at t = 1040/is for a lattice depth of Uo = 20Er 
and the initial 'coherent' state prepared by the a/6 displacement of the lattice. The echo amplitude is the ratio of the amplitude 
of the echo Gaussian envelope to 0.466, the amplitude of the original oscillations. A pu i se is the maximum lattice displacement, 
Wp U i Be is the temporal width of the pulse with T12 = 188/is at 20Er lattice depth, and P n >3 is the number of atoms in 
higher bands. The error bars for A pu i se and W pu i se indicate the uncertainties in our determination of the optimal pulse 
parameters, while those for the echo amplitudes are the direct measurement errors. The last column is the calculated P n >3 for 
the experimental pulse parameters. 



Due to the existence of original coherence, there is a re- 
phasing of oscillations and we observe an echo of the orig- 
inal oscillations centered at 2i = 2080/is. The echo am- 
plitude is extracted by subtracting the oscillations caused 
by the echo pulse itself. These oscillations are measured 
with each pulse applied at a time when there is no coher- 
ence left in the system and therefore we do not expect 
any revival of oscillations. In this experiment, each pulse 
is applied at t a = 4ms and the measured oscillations are 
then subtracted from the curves in Fig.[7Ja). Figure [T^b) 
shows the resultant echo curve after each pulse. 

The echo amplitude in the case of a single-step pulse at 
the earliest time we can measure, t = 2080/is, is w 0.19 
times the original oscillation amplitude, roughly an order 
of magnitude larger than the one observed in [20] . The 
main reason for this difference is that our lattice is an es- 
sentially conservative potential and the photon scatter- 
ing rate is negligible during the experimental run time 
(r sc = \QHz). In [2(|, the lattice is dissipative (cool- 
ing in the lattice) and atoms scatter a large number of 
photons (r sc = 2.5MHz). The lower echo amplitude in 
[2(j | is mainly due to photon scattering. By the time we 
measure the echo amplitude, i.e., at t = 2080/is ~ 11 os- 
cillation periods, the atoms have scattered 0.03 photons 
while in [201 ] . by the time the echo amplitude is measured, 
i.e., at t = 64/is ~ 12 oscillation periods, the atoms have 
scattered 160 photons. 

Some other principal differences in the two experi- 
ments are as follows. Our experiment is performed in 
relatively shallow lattices with lattice depths ranging 
18 — 20Er, while the lattice depths in 20J are in the 
range of 350 — 830-Er, where the coupling probability 
approaches the harmonic potential limit. We prepare 
atoms in the lowest band and measure the oscillation pe- 
riod between the lowest two bands while in |20j many 
bands are initially thermally populated. We measure the 
coherence between two energy levels and anharmonicity 
does not play a role in the decay of oscillations. In the 
Hannover experiment (20| many different anharmonic en- 
ergy levels are coupled through lattice displacements and 
the observed signal is a measure of mean atom position 
as the coherences between all of these levels evolve. In 
addition, the two experiments measure different things. 



In [2(| the authors measure the expectation value of the 
center of mass motion of the atoms, (x). We measure all 
state populations and can use this to perform complete 
reconstruction of the quantum state [29j]). 

Table [TT] gives a summary of pulse parameters, echo 
amplitudes and number of atoms lost from the lattice 
due to optimized pulses. In this case, each pulse is op- 
timized with respect to the amplitude of the echo. The 
envelope of the echo signal is fit to a Gaussian and the 
ratio to the original oscillation amplitude of 0.466 is pre- 
sented in table HU The square and Gaussian pulses result 
in larger echo amplitudes than the single-step pulse. The 
larger absolute echo amplitude after the Gaussian pulse 
is due to lower loss of atoms from the lattice than when 
the other pulse shapes are used. A normalization of the 
echo amplitudes with respect to the number of atoms left 
in the lattice (sum of the population of the two bands 
bound in the lattice) shows that the square pulse gives 
the largest relative echo amplitude. More atoms are lost 
from the lattice when we apply the echo pulse than when 
we prepare atoms in the lowest band and apply the opti- 
mized pulses. This is simply because at the time the echo 
pulse is applied there is a good fraction of atoms in the 
second band and the coupling from this band to higher 
bands is larger than the coupling from the lowest band 
to higher bands. For comparison to calculation, the last 
column in Table [TT] shows the calculated values for the 
coupling to higher bands. The calculated values show a 
similar trend to the experimental values in loss of atoms 
due to pulses. 

The amplitude and temporal widths of the experimen- 
tally optimized echo pulses are listed in Table IIIII and 
compared to the pulse parameters we calculate for opti- 
mum coupling in a lattice depth of 20-Er. The experi- 
mental echo pulse amplitudes are in agreement with cal- 
culation within the experimental uncertainties, but the 
agreement is not as good as in the case of the coupling 
experiment we discussed in section III Al (see Table UJ. 
The largest discrepancy is in the temporal width of the 
square pulse. The main difference between the two ex- 
periments, the coupling and echo experiments is in the 
initial conditions. In the pulse-echo experiment, at the 
time the pulse is applied there are atoms occupying both 
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Optimized Echo pulse Optimized Coupling pulse 



Pulse 


Apulse & Wp U l se 


P11 


P22 


P12 


Ap U lse & Wpulse 


P11 


P22 


Pa 


Single-step 


[0.29a, -] 


0.201 


0.052 


0.356 


[0.24a, -] 


0.334 


0.018 


0.388 


Square 


[0.167a, 0.40Ti 2 ] 


0.116 


0.095 


0.415 


[0.160a, O.32T12] 


0.234 


0.005 


0.470 


Gaussian 


[0.167a, O.3IT12] 


0.298 


0.020 


0.451 


[0.197a, O.25IT12] 


0.250 


0.013 


0.462 



Table III: Calculated coupling probabilities after single-step, square and Gaussian pulses for the pulse parameters used in the 
echo experiment (optimized echo pulse) and for the pulse parameters that give maximum Pi 2 (optimized coupling pulse) in 
a lattice depth of 20i?R. Pn and P22 are calculated considering all the population in the lowest or the next higher band. 
Pi 2 = P21 is the coupling probability between the two bands. 



lowest bands whereas in the case of coupling experiment 
initially all atoms are in the first band. An ideal echo 
pulse is the one that inverts populations of the lowest 
two bands. The optimum pulses discussed in section Hi Al 
result in maximum coupling to the second band, but they 
also leave a good fraction of atoms in the first band. 

Table IIIII lists the calculated number of atoms in each 
band after the single-step, square and Gaussian pulses 
for experimental optimized pulse parameters used in the 
echo experiment and compares those to the number of 
atoms in each band if the pulse parameters for optimum 
coupling P12 are used. P\\ and P22 are calculated cou- 
pling probabilities considering all atoms in band n = 1 
and band n = 2 respectively. P\i = P%\ is the coupling 
probability between the two bands. In the case of single- 
step and square pulses, with some compromise on the 
coupling, the pulse parameters used in the echo experi- 
ment leave fewer atoms in the lowest band and achieve 
a better population inversion and a larger echo ampli- 
tude. Figure [5] is a plot of the ratio P12/P11 to illustrate 
this fact. As can be seen in the figure, this argument 
does not hold for the Gaussian pulse. We have no in- 
tuitive explanation for this interesting observation. Fig- 
ure [5] also predicts when starting with atoms in the low- 
est band, after the experimentally optimized echo square 



pulse 78% of the remaining atoms will be in the second 
band (Pi 2 /(Pn + P12) = 0.78). Experimentally, we find 
this number to be f=a 0.63 after the optimized echo square 
pulse. 

We also measure how echo amplitude changes as a 
function of time. Figure [5] shows the measured echo 
amplitudes centered at 2tg after applying the optimized 
square pulse (for pulse parameters see Table [TTJ) at time 
U). Initially, the echo amplitude decays with an expo- 
nential decay constant of (720 ± 70) fis, then it stays ap- 
proximately constant over lms followed by a final decay. 
Further experiments are planned to understand this be- 
havior, which we believe provides information about the 
temporal correlation properties of the instantaneous well 
depth experienced by the atoms as they drift through a 
spatially inhomogeneous lattice. In the ID lattice, atoms 
are free to move in the transverse plane with an aver- 
age velocity, v rms cx s/T, where T is the temperature 
of atoms after molasses cooling. For example, a tem- 
perature of « 10/iK, a typical temperature after mo- 
lasses cooling, corresponds to v rms ~ 3.16cm/ s. An 
order of magnitude calculation for the time it takes for 
an atom with this transverse velocity to accumulate a tt- 
phase shift across a distance of 1mm, the rms radius of 
the lattice beams gives a decay time of ~ 1.7ms. 




Figure 8: (Color online) Bar plot of calculated ratio P12/P11 
for the pulse parameters that result in maximum P12 (opti- 
mized coupling pulse, checkered pattern bars) and the pulse 
parameters that result in maximum echo amplitude (opti- 
mized echo pulse, solid black bars). The blue bars show the 
experimental results. 
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Figure 9: (Color online) Echo amplitude as a function of time 
delay. The pulse used for the echo is the optimized experi- 
mental square pulse for the lattice depth of Uo ~ 20Er, and 
is centered at time t a ; the echo is observed at 2t . 
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(a) (b) 




q(7i/a) xnla 

Figure 10: (Color online) (a) Band structure of the periodic 
potential with Uo = ISEr, showing the first four bands. The 
dashed line indicates the depth of lattice. (b) ip ni g=o(xiT/a), 
the Bloch functions for quasi-momentum q — in the first 
four bands as a function of position. Lattice spacing is n with 
the central well extending from —n/2 to tt/2. The solid black 
curve is the lattice potential; the vertical dotted lines specify 
the center of wells and the horizontal dotted lines are drown 
for the illustration of wavefunctions. The Bloch functions 
have been rescaled and displaced vertically for clarity. 

III. MODELING 

In modeling our experiment, we start by calculating 
the eigenstates of the Hamiltonian H = P 2 /2m + U(x), 
where U(x) = U cos 2 (ki J x) is the periodic potential. 

h 2 k 2 

U = sEr is the lattice depth, where En = is 
the recoil energy, and k^ is the lattice vector. The 
eigenstates of the periodic potential are Bloch functions, 
^ n.q — e lqx u n-q (x), where q is the quasi-momentum and 
n is the band index 35]. Figure [TUT a) is the band struc- 
ture for a lattice depth of 18Er, a typical depth in our 
experiment. As an example of eigenstates, in Fig. [TWb) 
we plot the q — Bloch functions of the first four bands. 

Since in the experiment we prepare atoms in the low- 
est band of the lattice from a thermal cloud, all quasi- 
momenta of the lowest band are populated incoherently 
and with the same probability. This mixed state can be 
decomposed as an incoherent mixture of Wannier func- 
tions localized in each well. In one dimension, the Wan- 
nier function centered at x = is defined as [36[ 

w n ,l(x) = f ' dqip n , q (x)e~ tqla (2) 

J — ir / a 

Wannier functions are exponentially localized to each 
lattice site and they satisfy the relation w n .i(x) = 
Wn,i{x — a). In a tilted lattice, the eigenstates of the 
Hamiltonian H = P 2 /2m + V(x) + Fx are the metastable 
Wannier-Stark states, which form the discrete Wannier- 
Stark ladder [3(|. The Wannier-Stark states for the low- 
est lying states are localized to each well for small values 
of F. The lifetimes of the Wannier-Stark states are de- 
termined from the band-dependent Landau-Zener tran- 
sition rates ^lz = ^Bexp(—TT 2 aEQ/h 2 gn), where 
vb = Fa/h is the Bloch frequency, a is the lattice spac- 



ing, Eq is the band-gap, g is the acceleration due to grav- 
ity, and n is the band index for a transition from band 
n to band n + 1. For the lowest two states in our lattice 
depth of 18Er and the tilt per lattice site of 2.86-Er, the 
Landau-Zener transition rate for the n = 1 — * n = 2 
transition is ~ 4 x 10~ 7 Hz and for n = 2 — > n = 3 it is 
QHz. This results in lifetimes of 2.5 x 10 6 sec and 166ms 
respectively for these transitions, which are larger than 
the time scale of our experiment. For the n — 3 — > n = 4 
transition the rate is 1.2 x l0 3 Hz which corresponds to a 
lifetime of 830/is. Since both the Bloch oscillation period 
(530/is), and the Landau-Zener lifetimes are much longer 
than the temporal widths of the pulses used here, we be- 
lieve that the tilt of the lattice may be neglected. This 
is borne out by the agreement of our model with experi- 
ment, and has more recently been confirmed by compar- 
ison with more complete calculations including gravity 
(to be published elsewhere). 

A periodic Hamiltonian H = P 2 /2m + U cos 2 (kLx) 
conserves quasi-momentum. In addition, an abrupt 
displacement of the lattice does not change quasi- 
momentum. We therefore operate on each Bloch func- 
tion in the first band with a pulse O, then calculate 
the square of the overlap of the transformed Bloch func- 
tion with the Bloch function ^ n ^ q {x) of the same quasi- 
momentum in higher bands. Since we model the initial 
state as an incoherent sum over all quasi-momenta, we 
calculate the average transfer probability over all values 
of quasi-momentum, 




where 

\c ln \ 2 q = K*n,,(aOI<5|*i, 9 (aO>| a (3b) 

is the square of the matrix element of the operator O for 
each quasi-momentum. For example, to find the coupling 
due to a single-step pulse consisting of a lattice displace- 
ment, we calculate |ci„|^ = \( 1 ^ n ^ q (x)\ , ii, q (x — Acc))| 2 
for each quasi-momentum and integrate over all quasi- 
momenta to find the total coupling probability. In cal- 
culations, seven bands have been included and we have 
checked that including more bands (up to 20) does not 
change the result. 

Figure [TlTa) shows \ci2\ q , the coupling probability of 
the lowest two bands as a function of quasi-momentum 
for the lattice depth of U = 18Er. Figure fTlTb) is P\ n , 
coupling probability of the lowest band n = 1 with itself 
and with the second band n = 2 as a function of displace- 
ment. At this lattice depth, the coupling probability of 
the lowest two bands for the optimum displacement of 
Ax = 0.25a, with a the lattice spacing is P\2 = 0.39 
compared to the maximum coupling of 1/e « 0.37 for a 
harmonic potential. For comparison with the model, the 
measured values are shown on Fig. Illf b). The deviation 
of measured value for Pn from 1 and the measured value 
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Figure 11: (Color online) (a)|ci 2 |g = |(*n=2,«(a:)|*i, 9 (a: - 
Aa;)}| 2 as a function of quasi-momentum q for a single-step 
pulse (single lattice displacement) with the optimum displace- 
ment of Ax = 0.25a for the lattice depth of U a = 18Pr. (b) 
Pin = (1/2-Ka) f^i a / a dq\ci n \q as a function of displacement 
for bands n = 1 and n = 2 for the same lattice depth. For 
comparison, the measured values are shown on the same plot. 



for Pi 2 from at Ax /a — is due to transfer of atoms to 
the second band during the filtering stage (see section HT1 
and [1^1), which is done for the purpose of measuring 
populations of the two lowest bands. 

Each quasi-momentum contributes differently to the 
total coupling probability. Figure rT2] shows the coupling 
probability \c12\q for different single quasi-momentum 
components after single-step(a) and square(b) pulses as 
a function of displacement. 

The coupling probability changes as the lattice depth 
is changed. Figure fT3Fa) is the coupling probability Pi 2 
due to a single-step pulse as a function of displacement for 
different lattice depths. As the lattice depth is increased, 
the maximum coupling decreases and approaches the har- 
monic limit of 1/e at large lattice depths. To see the effect 
of anharmonicity of the lattice potential, we calculate the 
coupling probability for the harmonic potential. 

Figure [F5f b) is a plot of coupling probability Pi 2 as a 
function of displacement for different harmonic poten- 
tials. The frequency of the harmonic potentials have 
been chosen to match the average frequency of the low- 
est two bands of a particular lattice depth, i.e., we set 
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Figure 12: (Color online) |ci2 1, , the coupling probability due 
to (a) single-step pulse, and (b) square pulse(r = 0.35) for 
different quasi-momentum components as a function of dis- 
placement. 
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Figure 13: (Color online) (a) Optimum Pi 2 after a single-step 
pulse as a function of displacement for different lattice depths 
Uo = sEr. (b) The same overlap for harmonic potentials 
with cuho = IU12 , where ZU12 is the average frequency difference 
of the two lowest bands in a lattice depth of U a = sEr. Dis- 
placement in the case of harmonic potential is expresse d in 
units of lattice spacing a = ns 1 ^ 4 a, with a = ^Jh/mujho- 



LUho = S12 of the lattice. The maximum coupling in 
the harmonic case is 1/e and it is the same for differ- 
ent harmonic frequencies, but the magnitude of optimum 
displacement decreases as the harmonic frequency is in- 
creased. The coupling probability between the ground 
and excited states in a harmonic potential after displace- 
ment of the ground state can be calculated analytically. 
Using the momentum representation of wavefunctions for 
the ground and the first excited states, the coupling af- 
ter a single displacement of the ground state wavefunc- 
tion ]aJl\D(q)\0) = /+~ #(p)e-^to(p)dp = - 
Here, D(q) is the harmonic oscillator displacement oper- 
ator, £ = 221 i s a dimensionless displacement parameter; 



po = \j -^p is the width of the ground state wavefunc- 



tion in momentum space, (bo{p) 
and 4>i(p) = l , — ttt( — )e~ 



-e (2p °' 



r_2_y 



" (2^/4^- 

arc the momentum- 



space wavefunctions of the ground and first excited states 
of the harmonic oscillator. Then, the coupling probabil- 
ity is |(l|D(q r )|0) 1 2 = £ 2 e^ , which has a maximum of 
1/e when \q\ = = 2xo, where x = y / h/2muj is the 
width of the ground state wavefunction in real space. In 
fact, in the harmonic oscillator case any combination of 
displacements and time evolutions results in the maxi- 
mum coupling of 1/e. Since the optimum displacement 
xq is proportional to lu^ 1 ^ 2 as the harmonic frequency 
is increased the optimum displacement decreases. This 
general trend remains valid in the lattice. In the har- 
monic approximation of the periodic potential, the lat- 
tice potential U cx u> 2 , then the optimum displacement is 
proportional to P -1 / 4 . Therefore, as the potential depth 
is increased the optimum displacement decreases accord- 
ingly. 

To find the optimum parameters of a square pulse, 
we calculate the coupling probability with respect to the 
two parameters of the pulse, i.e., magnitude of the dis- 
placement and the time delay between the two displace- 
ments. The maximum overlap for a specific displacement 
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Figure 14: (Color online) (a) P12 after a square pulse as a 
function of displacement for different delays in a lattice of 
depth Uo = 18Er. t — delay /T12 is indicated in the legend 
for each curve. The optimum scaled delay between the two 
displacements for this lattice depth is r = 0.35 ( for U a = 
18E R , T 12 = 200/is, and optimum delay = 70/xs ). (b) P 12 
for the harmonic oscillator potential with uiho = £>i2 of the 
lattice. The delays in (b) are the same as in (a). Displacement 
in the case of harmonic potential is expressed in terms of 
lattice spacing a = ns 1 ^ 4 a, with a = ^h/musho- 



is found for different delays and the same calculation is 
repeated for different displacements to find the optimum 
displacement and the optimum delay. Figure I14f a) is a 
plot of Pi 2, the coupling probability of the lowest two 
bands as a function of displacement for different delays. 
For a lattice depth of U = 18Er wc find the maximum 
coupling probability of Pyi = 0.477 for the optimum dis- 
placement of Ax = 0.154a and the optimum delay of 
0.35Ti2, a 22% fractional improvement over the opti- 
mized single-step pulse (Table |T]). Again, we compare the 
lattice case with the harmonic potential. Figure E^b) is 
the plot of P12 due to square pulse as a function of dis- 
placement for a harmonic potential. In the case of a 
harmonic potential, the maximum coupling of 1/e stays 
the same for different optimum displacements and de- 
lays. Comparisons to harmonic potential in both cases 
of single-step and square pulses suggest that the anhar- 
monicity is the main reason for higher couplings in lower 
lattice depths. 

The interesting parameter in the square pulse is the 
delay, as this signifies how the phase difference between 
different states affects the coupling between the lowest 
two bands. The amount of phase accumulated during 
the time evolution in the lattice depends on the energy, 
and therefore the lattice depth. We calculate how the 
optimum delay changes as the lattice depth is changed. 
As the natural time scale of the system we choose Ti 2 , 
the average period of oscillations between the lowest 
two bands and present the scaled delay r as the ratio 
delay /Ti 2 . 

Figure[15]is a plot of the smallest optimum scaled delay 
r as a function of lattice depth. As the lattice depth is 
increased the scaled delay approaches 0.25Ti 2 . To have 
a better understanding of the optimum delay for deeper 
lattices, we calculate the optimum parameters of a square 
pulse for a harmonic potential with frequency lo^ = ZU 12 
of the lattice. Figure fTCT a) is a plot of the coupling P12 as 
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Figure 15: (Color online) The optimum smallest scaled delay 
of the square pulse, r = delay/Ti2 as a function of lattice 
depth. The optimum displacement is different for different 
lattice depths. 



a function of scaled delay r for three different fixed dis- 
placements of the harmonic potential. In the case of har- 
monic potential the maximum coupling stays the same 
for different displacements at optimum delays. In addi- 
tion there is no gain in coupling due to the square pulse 
over the single-step pulse. To see the effects of anhar- 
monicity alone in the case of lattice potential, we consider 
only one quasi-momentum component and observe how 
the coupling between the lowest two states changes as a 
function of delay. The plot in Fig. [TBT b) shows |ci2 |q = o 
as a function of r for the lattice depth of U — 18Pr. 
Similar to the harmonic potential, the optimum delay is 
different for different displacements, but unlike the har- 
monic case, the maximum coupling changes for different 
displacements. 

For a harmonic oscillator, one would intuitively expect 
the optimum delay to be T/2, half the period, since a 
displacement plus a delay of T/2 and a second displace- 
ment in the opposite direction is equivalent to a large 
single displacement. We draw a schematic of the phase- 



0.6 
0.5 
^0.4 
' 0.3 
0.2 
0.1 



(a) Harmonic Potential ■ Ax=o. 1 oa 

' .AX=0.16O 



(b) Lattice, q=0 



0.4 0.6 
T 



Figure 16: (Color online) (a) P12 after a square pulse as a 
function of scaled delay r for a harmonic potential with u^o = 
U12 in a lattice depth of U a = 18Er. Depending on the 
magnitude of displacement, the optimum delay changes but 
the maximum overlap P12 stays the same, (b) P12 after a 
square pulse as a function of scaled delay r for a single quasi- 
momentum component q = of the same lattice depth. P12, 
the maximum overlap is different for different displacements. 
Displacements in (b) are the same as in (a). 
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Figure 17: (Color online) (a) Phase-space schematics (not to 
scale) for the harmonic oscillator illustrating the square pulse 
sequence. The optimum delay t (8 — What) depends on the 
magnitude of displacement, (b) cos(8) = 1 — l/(2(Aa;/ri) 2 as 
a function of (Ax/n); ri is the phase-space radius of the first 
excited state of the harmonic oscillator. For the displacement 
of Aa; = 0.5ri the optimum delay is t — T/2 for cos(9) = — 1. 
As the displacement is increased the optimum delay becomes 
smaller; at Aa; = ri/\/2 (the blue dashed line), the delay of 
t — T/4 for 8 — tt/2 results in maximum coupling. 
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Figure 18: (Color online) Calculated Pi 2 as a function of 
lattice depth for optimized single-step, square and Gaussian 
pulses. The coupling probability increases as the lattice depth 
is decreased until the second band is no longer supported, at 
which part it drops off sharply. The inset shows Pi 2 for lattice 
depths less than 15Er. 



space representation of the harmonic oscillator number 
states in Fig. \T7\a,) to illustrate the effect of the dis- 
placement and time-evolution operators in coupling the 
ground state to the first excited state. In this represen- 
tation, each number state is represented as a circle of 
radius r n = y/n + 1/2 with the approximate width of 
1/2^/n. Since the number states are complete the cir- 
cles fill the phase-space. For sake of illustration the radii 
and the widths in Fig. [T7T a) are not drawn to scale. In 
this figure, the ground state corresponds to the central 
circle and the circular band represents the first excited 
state. A displacement corresponds to a translation along 
the x-axis and a delay or time evolution is a rotation in 
phase-space by an angle 6 = uit. As figure flTTa) shows, 
the optimum delay depends on the magnitude of displace- 
ment. In fact, for a displacement larger than T\j2 the op- 
timum delay is smaller than T/2. Figure fTTTb) is a plot 
of cos(0) = 1 — l/(Aa;/ri) 2 (calculated geometrically) as 
a function of displacement. The delay is t = T/2 only 
for a displacement of Aa; = 0.5ri. As the displacement is 
increased from r%/2 the optimum delay becomes smaller; 
for example at Ax = r\j\l2 the delay of t = T/4 for 
cos{6) — results in maximum coupling. 

Finally we compare the coupling probability for dif- 
ferent lattice depths due to optimized single-step, square 
and Gaussian pulses. We find that the square pulse yields 
a larger coupling probability than the single-step pulse 
for all lattice depths, with a fractional improvement as 
large as 30% for a lattice depth of U Q = 8Er. The cou- 
pling probability between the lowest two bands, P12, is 
shown as a function of lattice depth in Fig. [TH] for each 
of the three pulses. The effect of the Gaussian pulse is 
very similar to that of the square pulse except for lattice 
depths lower than 7 En, where the Gaussian pulse results 
in slightly larger coupling probability. For all pulses, the 
coupling increases as the lattice depth is decreased, until 
the second band becomes completely unbound, at which 



point it drops off sharply. The coupling probability for 
a single-step pulse increases from the harmonic potential 
limit of 1/e for large lattice depths to the maximum cou- 
pling of 0.51 for U = 5Tr. The square and Gaussian 
pulses perform even better, yielding couplings of 0.64 at 
6Er and 0.67 at 5Er, respectively. 



IV. CONCLUSION 

We have compared the effectiveness of three differ- 
ent pulses at coupling atoms in the lowest two bands of 
a far-detuned ID optical lattice. The three pulses are 
time-dependent lattice displacements following single- 
step, square and Gaussian profiles in time. We find 
experimentally that the square pulse results in larger 
coupling probability than the single-step and Gaussian 
pulses, a result which is in agreement with simulations. 
The maximum coupling probability observed experimen- 
tally, for a square pulse in a lattice of depth U — 18Er, 
is 0.33 ± 0.01; this is far higher than in any previous ex- 
periments, but lower than the theoretical maximum cou- 
pling of 0.477. We believe that this is due to lattice inho- 
mogeneity. We have also investigated optimizing pulses 
for inverting the populations of the lowest two bands in 
a pulse-echo experiment. We find that the square and 
Gaussian pulses result in larger echo amplitudes - 0.23 
and 0.24 times the amplitude of initial oscillations, re- 
spectively - than that arising from a single-step pulse, 
which is 0.19. The advantage of the Gaussian pulse arises 
because it leads to less loss from the lattice. For the op- 
timized echo square pulse, when starting with atoms in 
the lowest band 78% of the remaining atoms are pre- 
dicted to be transferred to the second band. Experimen- 
tally, 63% is our maximum observed value. We find that 
the agreement between experimental and calculated op- 
timized pulse parameters in the echo experiment is not as 
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good as the agreement in the case of the coupling exper- 
iment. The main difference between the two problems is 
that the echo pulse needs to efficiently invert populations 
from an arbitrary initial superposition of the two states. 
In addition, the pulses in the echo experiment have been 
optimized with respect to one parameter only, the echo 
amplitude for a particular initial coherence. A complete 
optimization of the pulses will require full quantum pro- 
cess tomography |29j ]. 

We have experimentally observed normalized echo am- 
plitudes of up to 0.24 in a lattice of depth 20Er, more 
than an order of magnitude larger than the echo am- 
plitudes observed previously for the motional states of 
atoms in optical lattices [20] ■ Our numerical calcula- 
tions show the coupling between the lowest two bands 
is more efficient in shallow lattices and should reach a 
maximum value of 0.67 for a Gaussian pulse when the 
lattice depth is 5 En, while a square pulse results in a 
maximum coupling of 0.64 for a lattice depth of 6Er. 
The square pulse outperforms the single-step pulse for 
all lattice depths and it results in larger couplings than 
the Gaussian pulse for lattice depths larger than 7Er. 
This is an example of coherent control, the optimized 
time delay introducing a relative phase which maintains 
constructive interference into the desired final state and 



destructive interference into the lossy states. These im- 
proved pulses have enabled us to extend pulse echo work 
to longer times, discovering a coherence plateau out to 
25 oscillation periods, a topic of ongoing study. 

Future work will include combining amplitude and 
phase modulation of the lattice potential to achieve bet- 
ter coupling and lower loss, investigation of adiabatic pas- 
sage approaches, as well as the application of optimal 
control theory to find improved pulse shapes and param- 
eters. The ability to efficiently couple quantized states in 
an optical lattice is crucial for state preparation as well 
as for the use of pulse echoes to preserve intra- well coher- 
ence, and similar pulses should find application in other 
systems described by tilted-washboard potentials. 
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